Anti-deSitter gravitational collapse 
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We describe a formalism for studying spherically symmetric collapse of the massless scalar field 
in any spacetime dimension, and for any value of the cosmological constant A. The formalism is 
used for numerical simulations of gravitational collapse in four spacetime dimensions with negative 
A. We observe critical behaviour at the onset of black hole formation, and find that the critical 
exponent is independent of A. 
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It is now well established that gravitational collapse in 
spherical symmetry exhibits a phase transition-like crit- 
ical behaviour, accompanied by self-similar behaviour of 
the matter field The basic formula determined 

numerically for the black hole radius near the threshold 
of black hole formation is p] . 



Rbh ~ (o - a*)~ 



(1) 



where a is an initial data parameter, and a* < a is its 
critical value. The critical value a = a* gives a naked 
singularity, and illustrates a violation of the cosmic cen- 
sorship conjecture. The exponent 7 is found to be "uni- 
versal" for a fixed matter type in the sense that it is 
independent of the functional form and parameter values 
of the initial matter profile. However it varies with the 
type of matter field. For perfect fluids for example, 7 
depends on the parameter(s) in the equation of stateQ. 

A wide class of matter types yield two similar general 
features near criticality: the black hole radius scaling law 
(without a mass gap), and discrete or continuous self- 
similarity of the matter fields. (These are called Type II 
transitions. There are also transitions with a mass gap, 
refered to as Type I, where the field is either periodic or 
static, rather than self-similar.) 

One of the questions that has not been studied to date 
is that of the dependence of the critical exponent on the 
cosmological constant A. We study this question in four 
spacetime dimensions. (The case of three dimensions is 
rather special in that there are no static black hole so- 
lutions unless the cosmological constant is negative; the 
constant may be scaled to unity without loss of general- 
ity so there is no need to study such dependence in this 
case d,@)- Our main result is that 7 is independent of 
A. 

A useful by-product of our analysis is confirmation of 
the utility of a new formalism for studying the collapse 
problem, in which the d— dimensional, spherically sym- 
metric Einstcin-scalar equations are rewritten as an effec- 
tive 2— dimensional dilaton gravity theory 0- All space- 
time and A information in this approach is stored in the 



resulting dilaton potential. This feature permits writing 
a "universal" numerical code in which spacetime dimen- 
sion and A appear as input parameters. In fact the for- 
malism is general enough to permit simulations of mass 
and potential terms for the scalar field, with minimal 
changes to the code. 

The reduced equations are written in double null co- 
ordinates, and a numerical method first used by Gold- 
wirth and PiranQ, and refined by Garfinkle|l(| is imple- 
mented on the resulting equations. The code is in prin- 
ciple capable of handling any spacetime dimension and 
cosmological constant value, although there are practical 
constraints. We describe below the main features of the 
formalism and numerical method. Further details appear 
in Ref. 0. 

Einstein gravity with cosmological constant and min- 
imally coupled scalar field in d spacetime dimensions is 
given by the action 
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(2) 



Spherical symmetry is imposed by writing the metric g av 



as 



ds td) = g a fsdx a dx (i + r 2 (x a )dtt( d - 2 ), 



(3) 



where dCl/4-2) is the metric on S d ~ 2 and a, (3 — 1, 2. 

A useful form for the reduced action with this form of 
the metric is obtained by defining I — (G^ d ') n ' 2 and 



8(n-l) 



la/3, 



(4) 
(5) 



where n = d — 2. Note that the <f> is proportional to 
the area of an n-sphere at fixed radius r. With these 
definitions the reduced action becomes 



s - JL 

S ~ 2G 



d 2 Xy 
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4>R(g) + vW(ct>,A) 
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- / d 2 xV=g H (n \cj>) \d x \ 2 (6) 

where G = 2nn/ {n— 1), the scalar field has been rescaled 
(x ~~ * x/vG^) in order to make it dimensionless, and 
the overall factor of the n— sphere volume has been 
dropped from the action. In addition we have defined 



n 



A) 



8(n - 1) 



8(n- 1) 



Now with the metric parametrized as 

ds 2 = —2lg(u,v)(f>'(u,v)dudv 
the field equations are 



g'<t>' 



2G( X 



/\2 



(7) 

(8) 

(9) 

(10) 
(11) 



(»)(</>) 

(ffW(^)' + (ff( B) Wx)' = 0, (12) 

where prime and dot denote the u and u derivatives re- 
spectively. 

The evolution equations may be put in a form more 
useful for numerical solution by defining the variable 



X + 



20X' 



(13) 



which replaces the scalar field x by /i. The evolution 
equations become 



-5/2 
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where 
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An J dv^{h-x? 
(g<t>'VW)dv, 
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2v^ Ju 



dv 



h<t>' 



(14) 
(15) 

(16) 
(17) 

(18) 



This is the final form of the equations used for numerical 
evolution. The spacetime parameters n and A appear 
only in the dilaton potential V^ 1 ' ©• 

The numerical scheme uses a v discretization to obtain 
a set of coupled ODEs: 



where % = 0, • • • , N specifies the v grid. Initial data for 
these two functions is prescribed on a constant v slice, 
from which g(u, v), g(u, v) are constructed. Evolution in 
u ('time') is performed using the 4th. order Runge-Kutta 
method. The general scheme is similar to that used in 
0, with some refinements used in [Tof . This procedure 
was also used for the 3— dimensional collapse calculations 
in 0. 

The initial scalar field configuration %(</>, u = 0) is most 
conveniently specified as a function of <fr rather than r. 
(Recall that <f> oc r n .) This together with the initial ar- 
rangement of the radial points <p(v,u — 0) fixes all other 
functions. We used the initial specification 0(0, v) = v. 

Most of the computations are for the initial scalar field 
of the Gaussian form 



Xg{u = 0, (f>) = a<j) exp 



(20) 



with attention restricted to variations of the amplitude 
a. However we also used the "cosh" initial data 



Xg(u = 0,(p) = acosh[b((f) - O )]~ 



(21) 



to study convergence of our code and to test universality 
with non-zero A. 

The initial values of the other functions are determined 
in terms of the above by computing the integrals for g n 
and g n using Simpson's rule. 

The boundary conditions at fixed u are 



= 0, <7fe=0, gk = 1. 



(22) 



where k is the index corresponding to the position of the 
origin (f> — 0. (In the algorithm used, all grid points < 
i < k— 1 correspond to ingoing rays that have reached the 
origin and are dropped from the grid). These conditions 
are equivalent to r(u,u) — 0, g\ r =o — g{u,u) — 1, and 
guarantee regularity of the metric at r = 0. For our initial 
data, 4>k and hence hk are initially zero, and therefore 
remain zero at the origin because of Eqn. (|15|) . 
At each u step, the function 



ah = g a0 d a 



ig' 



(23) 



is observed. Its vanishing signals the formation of an 
apparent horizon. For each run of the code with fixed 
amplitude a, this function is scanned from larger to 
smaller radial values after each Runge-Kutta iteration, 
and evolution is terminated if the value of this function 
reaches 10 -4 . (The code cannot continue much beyond 
this value.) The corresponding radial coordinate value is 
recorded as R a h- In the subcritical case, it is expected 
that all the radial grid points reach zero without detec- 
tion of an apparent horizon. This is the signal of pulse 
reflection. The results (a,R a h) are collated by seeking a 
relationship of the form 



h(u, v) — > hi{u), 



(u,v) -> (j)i{u). 



(19) 



R ah oc (a - a») 7 . 



(24) 
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The code was tested for grid sizes ranging from 2000 
to 20000 points, and with the u and v step sizes ranging 
from 10 -2 to 10~ 6 , for the two types of initial data used, 
as well as the vacuum case x = 0. These tests established 
that the code converges. 

For most of the runs we used <po = 1 and a = 0.3 for 
the Gaussian initial data, and b = 5.0, </>o = 0.5 for the 
cosh data. We varied <po in the Gaussian data from 0.1 
to 20 (which also varies pulse width in the r-coordinate) , 
to see if there was 7 dependence. This is important to 
study because A sets a scale in the problem. 

Our results for the two types of data (for n — 2) appear 
in Tables 1 and 2, which list the computed 7 values for 
the corresponding (negative) values of A. The trend is 
clear: 7 does not depend on A. The first value provides a 
crucial test of our formalism and code since it reproduces 
the known A = value of the exponent 0. 



A 7 

-0.001 0.370-0.375 

-5 0.37-0.38 

-10 0.37-0.39 

-20 0.36-0.38 



Table 1. Computed ranges of 7 for Gaussian data (n = 2). 



A 7 

^5 0.37-0.38 

-20 0.36-0.38 

-50 0.37-0.40 



Table 2. Computed ranges of 7 for cosh data (n = 2). 

The ranges of 7 in Tables 1 and 2 are arrived at by as- 
sessing our uncertainty in the determination of the crit- 
ical amplitudes a*: the a* values are determined to lie 
within certain domains, and the end points of these do- 
mains are used to determine the range of 7 values. For 
example, for A = — 20 and cosh data, we find a* lies in 
the range 6.00344 x 10~ 3 - 6.00360 x 10" 3 , which gives 
7 in the range indicated in Table 2. Our determination 
of the a* window gets coarser with increasing A, which 
leads to the larger error bars on the 7 values, as indicated 
in the tables. 

Figure 1 is the ln-ln graph of the apparent horizon 
radius vs initial data amplitude, for A — —5. This is rep- 
resentative of the type of results from which the critical 
exponent ranges in Tables 1 and 2 are deduced. 

Figure 2 shows plots of the apparent horizon function 
(|23*|l and the scalar field h. This is typical of the type 
of data used to extract horizon radius information at the 
supercritical onset of black hole formation. These graphs 
show impending black hole formation at <\> ~ 6.3 x 10~ 6 
(r ~ 5 x 10~ 3 ) for a supercritical amplitude with A = —5. 
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In(a-a-) 

FIG. 1: Logarithmic plot of apparent horizon radius R a h ver- 
sus initial scalar field amplitude (a — a*) for A = —5. The 
line is the least squares fit to the points giving 7 = 0.3738 
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FIG. 2: Typical graphs of the apparent horizon function (I2.'-!I 
and the scalar field h near black hole formation. 



The dip in the apparent horizon function (Fig. 2) os- 
cillates toward and away from the origin as it approaches 
the <fi— axis, and appears to be the cause of the small 
oscillations about the least squares fit line in Figure 1. 
Although these latter oscillations have been noted in ear- 
lier work, it appears that they are the manifestation of 
the dip oscillations in the apparent horizon function 1231) 
near criticality. This in turn is connected with the dis- 
crete self-similarity of the scalar field, which not surpris- 
ingly, is also manifested in the apparent horizon function. 
Figure 1 contains two complete oscillations, and may be 
used to obtain a rough estimate of the echoing period: 
it gives A ~ 3.5. This value is close to the A = value 
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A = 3.44 computed in earlier workQ, 0, Thus, both 
the numbers 7 and A associated with critical behaviour 
appear to be independent of A. 

It is known, and verified again here, that 7 is indepen- 
dent of initial data profiles and parameter values. We 
also find that for fixed A, critical exponents are indepen- 
dent of characteristic pulse width and starting location 
in relation to the A scale. 

It is worth noting that in our equations A appears only 
in the potential term JSJ, where for sufficiently small </>, 
the A term is insignificant. This provides some analyti- 
cal support of our results, since it suggests that the entire 
near critical part of the evolution, with typical <fr range 
— 10~ 5 (Fig. 2), is A independent. In four dimensions, 
n = d — 2 = 2, the first term in brackets in (jHJ is of order 
10 5 , which dominates the A term for all cases we consider. 
An extrapolation of this observation, beyond numerical 
reach, suggests that 7 is A independent because suffi- 
ciently close to criticality, the <j) term dominates the A 
term for any A (in an e — 6 sense) . 

It is interesting to contrast this with the minimally 
coupled scalar field of mass /1, which also has a scale. 
Here two types of critical behaviour are observed at the 
threshold of black hole formation, de pen ding on initial 
pulse width a in comparison with \i~ x 12 : for a >> /i -1 
there is a mass gap, whereas for a << /i -1 there is no 
mass gap and the exponent 7 ~ 0.378 is computed [l2|. 

There are some differences here worth emphasizing. 
The first is specific to our results, and the second is gen- 
eral: (i) We find no evidence of a mass gap for a range 
of values of pulse width and initial position in relation to 
the scale l/V— A. With no mass gap, the critical expo- 
nent is measured in the limit of small horizon values. In 
our simulations, the typical horizon size range at the on- 
set of black hole formation is 10~ 2 — 10~ 4 . By contrast, 
the smallest A scale the code can handle is 1/V50 ~ 0.14. 
Under these circumstances, it is not surprising that A has 



no effect, (ii) Unlike (1, A appears in the stress-energy 
tensor independent of the scalar field. Therefore it is not 
unreasonable that A and /j, give qualitatively different 
results. 

For postive A we are not able to obtain accurate results 
for 7 because of two competing effects in our procedure: 
since the grid is evolving, we observe an outward deSitter 
expansion of the grid, which confines the interesting fea- 
tures of the ingoing collapse to an ever shrinking region 
near the origin. Perhaps this case can better studied by 
the method of replacing lost grid points used in |Tcj , and 
using the subcritical method of computing 7 using the 
curvature scalar, introduced in [l^. Nevertheless, on the 
intuitive grounds mentioned, we expect that the sign of 
A will not change our results. 

It is known that critical exponents may be computed 
via linear perturbation analysis of the critical solution 
0,^3 ■ This has so far been done only for A = 0. For A ^ 
0, it is reasonable to expect that both the critical solution 
and its perturbation equation depend on A, which after 
all is a parameter in the equations of motion. It would be 
useful to see how A drops out of this type of calculation, 
yielding a A independent exponent. 

In summary, our numerical simulations of massless 
scalar field collapse in spherical symmetry in four dimen- 
sions, show that the critical exponent associated with the 
collapse is independent of (negative) A values. This result 
extends the scope of universality to include the cosmolog- 
ical constant, and suggests that including A with other 
matter types, such as the perfect fluid or the Yang-Mills 
field, will also not change the critical exponent. 
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